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Abstract 

The decay of the Z' boson into supersymmetric particles is studied. We investigate how these 
supersymmetric modes affect the current hmits from the Tevatron and project the expected sen- 



r~ I ! sitivities at the LHC. Employing three representative supersymmetric Z' models, namely, Eq, 

U{1)b-l, and the sequential model, we show that the current limits of the Z' mass from the Teva- 



tron could be reduced substantially due to the weakening of the branching ratio into leptonic pairs. 

The mass reach for the Eq Z' bosons is about 1.3 — 1.5 TeV at the LHC-7 (1 fb~^), about 2.5 — 2.6 
(N 
> ■ TeV at the LHC-10 (10 fb"^), and about 4.2 - 4.3 TeV at the LHC-14 (100 Ib'^). A similar mass 

^^ ' reach for the U{1)b~l Z' is also obtained. We also examine the potential of identifying various 

supersymmetric decay modes of the Z' boson because it may play a crucial role in the detailed 

^^ I dynamics of supersymmetry breaking. 



X 



I. INTRODUCTION 

Existence of extra neutral gauge bosons has been predicted in many extensions of the 
standard model (SM) [l|. String- inspired models and grand- unification (GUT) models usu- 
ally contain a number of extra U{1) symmetries, beyond the hypercharge f/(l)y of the SM. 
The exceptional group Eq is one of the famous examples of this type P]. These extra U{l)'s 
are broken at some intermediate energy scales between the GUT and the electroweak scales. 
Phenomenologically, the most interesting option is the breaking of these U{l)^s at around 
TeV scales, giving rise to extra neutral gauge bosons observable at the Tevatron and the 
Large Hadron Collider (LHC). Recent developments in model buildings also result in new 
models that contain extra gauge bosons. For example, little Higgs models J3| with additional 
gauge groups predict a number of new gauge bosons; the SM gauge bosons propagating in 
the extra dimensions after compactification can give rise to Kaluza-Klein towers of gauge 
bosons |j]; Stueckelberg Z' model connecting a hidden sector to the visible sector in the 
context of dark matter [5|, just to name a few. We will denote these neutral extra gauge 
bosons generically by Z' for the general discussion in this introduction section. 

Collider experiments such as CDF 6| and D0 [7|] at the Tevatron have been searching for 
the neutral extra gauge bosons Z', mainly through its leptonic decay modes. The leptonic 
mode is a very clean channel to probe for Z' since it may give rise to a discernible peak 
above the Drell-Yan background right at the Z' mass, provided that the size of the coupling 
strength to SM quarks and leptons is not too small. Currently, the best limit comes from 
the negative search at the Tevatron. The lower mass bound on Z' is about 800 — 900 GeV 
for a number of Z' bosons of the Eq type and a stronger bound of almost 1050 GeV for the 
sequential Z', which has exactly the same coupling strength and chiral couplings as the SM 
Z boson that can be served as a bench mark. The LHC with just an integrated luminosity of 
about 40 pb~^ has already set limits on Z's 8| almost as good as those from the Tevatron. 
With more luminosity accumulated in the current LHC run the limit on Z' will improve 
substantially in the near future. 

Most of previous studies on Z' bosons focused on the decays into SM fermions, and 
the corresponding limits were obtained based on the decay into leptons. This scenario is 
not necessarily a must, but just for simplicity and fewer choices of parameters. Indeed, 
when the mass of Z' is more than a TeV or even larger it has chances of decaying into 



other exotic particles that must be included in the model for various theoretical reasons. 
For instance, in GUT models or in little-Higgs models there are other fermions needed 
to cancel the anomalies, the Z' could decay into these exotic fermions if their masses are 
not too heavy. Another example is the minimal supersymmetric standard model (MSSM) 
with electroweak-scale SUSY partners and a Z' which could also decay into sfermions and 
Higgsinos. In a recently proposed Z'-mediated SUSY-breaking model 9|, supersymmetry 
breaking in the hidden sector is communicated by a Z' boson to the visible sector. The 
low-energy spectrum includes a Z' boson of a few TeV and light gauginos, such that the Z' 
can also decay into SUSY particles other than the SM fermions. 

In this work, we consider a scenario of Z' boson, which arised from f/(l) symmetry 
breaking at around TeV scale, in the context of weak-scale supersymmetry, in which all 
SUSY partners are relatively light (a few hundred GeV) except for the squarks (may of 
order 0(1) TeV) . Such a Z' may come from breaking of one of the ?7(l)'s in iJg (2[, U{X)b-l 



10|, or f/(l)_B 111 ^^"^^ Oiice we specify the f/(l) charges for the matter superfields and Higgs 
superfields, the couplings of Z' to MSSM particles are determined. We study the decays of 
Z' and its production at the Tevatron and the LHC. The decays of Z' will be modified when 
the SUSY particle masses are only of order a few hundred GeV, which will then affect the 
leptonic branching ratio of the Z' . We investigate how much the limits from the Tevatron 
will be affected, because of the reduction in the leptonic branching ratio. We also study how 
much the sensitivity at the LHC will be reduced when the SUSY decay modes are open for 
the Z' boson. Finally, we study the prospect of using the SUSY decay modes of the Z' to 
search for the Z' boson itself and investigate its properties. 



We note that some related works had appeared in literature, for example in Refs. 12Hl4|. 
However substantial improvements over these previous works have been made in this work. 
These include 



1. Ref. 12| focused on how the presence of supersymmetric and other exotic particles of 
E^ models in the Z' decay can affect the Z' boson discovery at the Tevatron and the 
LHC. In our work, we used the most current updated limit on cr(Z') x Bi^Z' — )■ fi~^fi~) 
to put limits on the mass of various Z' bosons. We have illustrated the case of decaying 
into SM particles only and the case of including both SUSY and SM particles. We 
have shown that the Z' mass limits have to be relaxed by 20 — 30 GeV if including 



SUSY particles in the decay. Ref. 12| was written in 2004 and certainly our paper 
used the newest 2011 data. For the LHC sensitivity we worked out the more realistic 
energy- luminosity combinations (7 TeV, 10 TeV, 14 TeV). Refs. 
about the options of 7 and 10 TeV at their time. 
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13| did not know 



2. Ref. 13[ focused on how the discovery potential of sleptons can be improved via the 
decay of Z' into a slepton pair. They also studied various lightest supersymmetric 
particle (LSP) scenarios and used distributions to determine the masses of sleptons, 
gauginos, and the Z' boson. In our work, in addition to the slepton-pair, we also 
studied the Z' decays into a chargino pair and into a neutralino pair. We have shown 
clearly that the presence of Z' is visible in the transverse-mass spectrum. One can 
therefore measure the transverse-mass spectrum and determine if there is a Z' boson. 
This spectrum can also be utilized to estimate the mass differences, and couplings of Z' 
to sleptons, neutralinos, and charginos. This can help us to understand the underlying 
supersymmetry breaking mechanism. 



3. In this work, we study E^, models, f/(l' 
studied only the E^ models and Ref 



R /-,, and the sequential Z' model, while Ref. [12 1 
13| studied only the U{X)b~xL model. 



The organization of the paper is as follows. In the next section, we write down the inter- 
actions and briefly describe a few Z' models and their extensions to include supersymmetry. 
In Sec. Ill, we calculate the branching ratios of the Z' boson in various models. In Sec. IV, 
we show the shift of the limits for the masses of the Z' in various models due to opening 
of supersymmetric particles. We estimate the Sa discovery reach at the LHC, including the 
SUSY decay modes in Sec. V. We further discuss in Sec. VI the SUSY decay modes of the 
Z' boson. We conclude in Sec. VII. Feynman rules that are related to the Z' are collected 
in the appendix. 



II. Z' INTERACTIONS AND REPRESENTATIVE MODELS 

A. Z' Interactions 

Following the notation of Ref . [l|, |l5| , the Lagrangian describing the neutral current gauge 
interactions of the standard electroweak SU{2)l x U{1)y and extra f/(l)'s is given by 

n 

- £nc = eJ^^A, + Y, g.JliZl , (1) 

a=l 

where Z^ is the SM Z boson and Z° with a > 2 are the extra Z bosons in the weak-eigenstate 
basis. For the present work we only consider one extra Z^ mixing with the SM Z^ boson. 
Thus the second term of the Lagrangian in Eq. ([T]) can be rewritten as 



^iM 



Y,i^n'i9iPL + giPR)i^f 



+ 92ZI 



^zozo = giZi 

L / J L / J 

(2) 



Y,i^friQ'f,PL + Q'f,PR)i'f 



where for both quarks and leptons 

9l,r — -^3L ~ ^wQf , (3) 

and Q'f r are the chiral charges of fermion / to Z2 and Pl,r = (1 =F 75)/2. Here T^^ and 
Qf are, respectively, the third component of the weak isospin and the electric charge of the 
fermion /. The chiral charges of various Z' models are listed in Tables HI and ITTl The overall 
coupling constant gi in Eq. ([1]) is the SM couphng g/ cosO^, while in grand unified theories 
(GUT) (72 is related to gi by 

^=(^x.Ay^'^0.62A^/^ (4) 

where x^ = sin^ 6„ and 6„ is the weak mixing angle. The factor A depends on the symmetry 
breaking pattern and the fermion sector of the theory, which is usually of order unity. 

The mixing of the weak eigenstates Z^ and Z2 to form mass eigenstates Z and Z' are 
parametrized by a mixing angle 6: 




cos 6 sin 6 
- sin 6 cos 6 




(5) 



The mass of Z is Mz = 91.19 GeV. 



After substituting the interactions of the mass eigenstates Z and Z' with fermions are 



/ 



Z,i,a\l{P^ + r{Pn)i^^ + Z'^a\l{P^ + r{Pn)^j 



(6) 



where 



'{=9{ + |#Q',,, ri=g<„+'f^eQ)^, (7) 

li = 'l^i-09{. ri='fQ'i-e,i. (8) 

Here the subscript "s" denotes the observed SM Z boson and "n" denoted the new heavy 
gauge boson Z' . We have used the vahd approximation cos^ ~ 1 and sin^ f^ Q. In the 
following, we ignore the mixing [Q = 0) such that the precision measurements for the SM Z 
boson are not affected, unless stated otherwise. 



B. Including Super symmetry 

The superpotential W involving the matter and Higgs superfields in a t/(l)' extended 
MSSM can be written as 

W = ea, [y:t,Q]HtU^ - vt.QW'dD- - y\,^H',Et + hsSH^H^] , (9) 

where ei2 = — £21 = 1? hj are family indices, and y^ and y'^ represent the Yukawa matrices 
for the up-type and down-type quarks respectively. Here Q, L, If^, D'^, E'^, H^, and H^ denote 
the MSSM superfields for the quark doublet, lepton doublet, up-type quark singlet, down- 
type quark singlet, lepton singlet, up- type Higgs doublet, and down-type Higgs doublet 
respectively, and the S is the singlet superfield. Note that we have assumed other exotic 
fermions are very heavy. The f/(l)' charges of the fields Hu, H^, and S are related by 
Q'hu + Q'h + Q's ~ ^^'^'^ ^^^^ SH^Hd is the only term allowed by the U{\)' symmetry 
beyond the MSSM. Once the singlet scalar field S develops a VEV, it generates an effective 
/i parameter: /icg = hs{S). The case is very similar to NMSSM, except we do not have the 
cubic term 5*^. The singlet field will give rise to a singlet scalar boson and a singlino, which 
will mix with other particles in the Higgs sector and neutralino sector, respectively. In this 
work, we are contented with the assumption that the singlet scalar field and the singlino 
are heavy enough that it is out of reach at the LHC and the mixing effects are negligible. 
Detailed phenomenological studies involving the singlet field will be presented in a future 
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work. Furthermore, we also take the superpartner, dubbed as Z'-ino, of the Z' boson to be 
heavy. Phenomenology involving the singlet scalar boson, singlino, and the Z'-ino of various 
singlet-extended MSSM can be found in Refs. [l6|. 

Below the TeV scale the particle content is the same as the MSSM plus a Z' boson. Thus, 
the superpotential and the soft breaking terms are the same as in MSSM. Extra couplings of 
the Z' boson with the MSSM particles are coming from the gauge interactions of the extra 
f/(l) and the corresponding supersymmetric vertices of Yukawa interactions. 

The gauge interactions involving the fermionic and scalar components, denoted generi- 
cally by ip and respectively, of each superfield are 

C = ^ t^D, ^ + {D^y {D^(t>) , (10) 

where the covariant derivative is given by 



D^ = d, + leQA^ + ^^(^+^,t + ^"^m") + ^^i(^3L - Qx^)Z^ + ig^Z'^Q' . (11) 

Here e is the electromagnetic coupling constant, Q is the electric charge, r^ are the rising and 
lowering operators on SU(2)l doublets and Q' is the chiral charges of the f/(l) associated 
with the Z' boson. The interactions of Z' with all MSSM fields go through Eqs. flTUl) and 

m- 

Details and conventions are given in the appendix and the Feynman rules that involve 
the Z' boson are listed there as well. 



C. Representative models 

1. The Sequential model 

The sequential Z'gj^j model is a reference model of extra Z bosons. It has exactly the 
same chiral charges as the SM Z boson but at a larger mass. The gauge coupling constant 
is also taken to be the same as the SM one, i.e., fi'2 = fi'i = g/ cos9^. Note, however, that 
when SUSY modes are open, the Zgj^ can also decay into sfermions, neutralinos, charginos, 
and Higgs bosons. In general, experimental constraints on the sequential model are the 
strongest, because the other Eq models, for example, have smaller gauge coupling constant, 
as in Eq. (jl). 



2. The Eq models 

Two most studied f/(l) subgroups in the symmetry breaking chain of E^ occur in 

Ee -^ 50(10) X U{1)^ , 50(10) -^ SU{5) x U{1)^ . 

In Eq each family of the left-handed fermions is promoted to a fundamental 27-plet, which 
decomposes under Eq — )■ 5*0(10) — )■ SU{5) as 

27 -> 16 + 10 + 1 ^ (10 + 5* + 1) + (5 + 5*) + 1 . 

Each 27 contains the SM fermions, two additional singlets u'^ (conjugate of the right-handed 
neutrino) and S, a. D and D'^ pair {D is the exotic color-triplet quark with charge —1/3 
and D'^ is the conjugate), and a pair of color-singlet SU(2)-doublet exotics H^ and Ha with 
hypercharge Yh^h^ = ±1/2. In the supersymmetric version of Eq, the scalar components of 
one Hud pair can be used as the two Higgs doublets H^d of the MSSM. The chiral charges 
f/(l)^ and f/(l)x for each member of the 27 are listed in the third and fourth columns in 
Table [H In general, the two U{1)^ and U{1)^ can mix to form 

Q'i^Ee) = cos 9 EeQ'^ + sm 9 EsQ'^i, , (12) 

where < 6'^g < vr is the mixing angle. A commonly studied model is the Z' model with 

Q'n = xHq'. - \IIq'^ ' (13) 



which has 9eq = tt — tan ^ a/5/3 ~ O.Tlvr. There are also the inert model with Q'j = 
-Q'{9e^ = tan~^ ^/3/5 ~ 0.2l7r), the neutral N model with 9Ee = tan~^ v^ ~ 0.427r, and 
the secluded sector model with 9Ee = tan~^ a/15/9 ~ 0.137r. The chiral charges for each 
member of the 27 are also listed in the last four columns in Tabled for these four variations 
of Z' models within Eq. Here we take the assumption that all the exotic particles, other 
than the particle contents of the MSSM, are very heavy and well beyond the reaches of all 
current and planned colliders. 

3. B — L models 

The extra U{1) symmetry here is the U{1)b^l with Q'f = ^{B — L)^ , where B and 
L are the baryon and lepton numbers, respectively. The theory is vector-like, thus the 
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TABLE I. The chiral charges of the left-handed fermions for various Z' bosons arised in i^g Ul- 
Note that Q'r = —Q'if^) since aU the right-handed SM fermions are necessarily converted into 
left-handed charge-conjugated fermions in order to put them into the irreducible representation of 
27 of Eq. 
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chiral charges for the left- and right-handed fermions are therefore the same, Q'^ = Q'f = 
^(B — L)f as they are summarized in Table [ITl 

TABLE IL The chiral charges of the left-handed fermions for the Z'^_^ boson. Note that Q'r = 
—Q'if^) since charge-conjugated fields are used here. 

Q'b-l = {B- L)/2 

Q = {u,d) 1/6 

u"" -1/6 

d^ -1/6 

L = {u,e-) -1/2 

e^ 1/2 



III. DECAYS OF Z' 

In general, the Z' boson can decay into all SM fermion, squark, slepton, sneutrino, neu- 
tralino, chargino, and Higgs-boson pairs. The decays into neutralino and chargino pairs go 
via the couplings to Higgsinos. The Higgs-boson pairs include Ah'^, AH^, and H^H^ . All 
decay modes are subjected to kinematic threshold. We have chosen typical input parameters 



for the supersymmetric particles: [j |i 



Set (A) : tan /3 = 5 , mg = 0(a few TeV) , Mi = 200 GeV , Ma = 500 GeV , 

Ml = 100 GeV , Ms = 200 GeV , /i = 150 GeV , A = 100 GeV . (14) 

With this set of choices the Z' boson will not decay into squark pairs, but can decay into 
all other sfermion pairs, neutralino and chargino pairs, and Higgs-boson pairs. 

We show in Figs. [1] to [8] the decay branching ratios of various Z' bosons. In these figures, 
we only show e~^e~ - one of the charged lepton modes, which is the most direct discovery 
mode of the Z' boson, e^e^ + e/je^j - one of the slepton modes, UeL^lh ^ ^^^ of ^^^ sneutrino 
modes, but the sum of neutralino and chargino pairs. The other charged-lepton modes and 
slepton modes are approximately the same as the corresponding one shown in the figures. 
The top-quark pair and the other Higgs-boson pairs are also shown. Note that the Z'^_]^ 
does not couple to the Higgs fields. 

We also select another set of SUSY parameters that we call set (B), typically it is a large 
tan/3 case. 

Set(B) : tan /3 = 40 , m^ = 0(a few TeV) , M^ = 300 GeV , Ma = 500 GeV , 

Ml = 200 GeV , Ms = 400 GeV , /i = -300 GeV , A = 300GeV . (15) 

We shall show the results for the limits from the Tevatron and the sensitivities at the LHC 
with the two choices of set (A) and set (B). 



^ There are additional _D-term contributions to the sfermion and Higgs boson masses from the breaking of 
the U{iy symmetry [17|, which further depend on the details of the VEV of additional Higgs fields that 



break the [/(I)' symmetry. Instead, we used the physical masses as inputs to our study. 
^ We choose the squark mass to be heavy enough such that the decay into squark pair is not open. We 

have set the soft parameters A/? = Mf^^ = (200 GeV)^ such that the physical slepton masses are about 

205 and 204 GeV for the left-handed and right-handed sleptons, but about 190 GcV for the sneutrino. 

We have ignored the additional D-term contribution from the 1/(1)' . 
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FIG. 1. Decay branching ratios for some decay modes of the sequential Z' modeh 
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FIG. 2. Same as Fig. [D but for Z' model. 
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FIG. 3. Same as Fig. [T]but for Z'. model. 
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FIG. 4. Same as Fig. [U but for Z' model. 
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FIG. 5. Same as Fig. [T] but for Z'j model. 
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FIG. 6. Same as Fig. [T]but for Z'^ model. 
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FIG. 7. Same as Fig. [T] but for Zg^^ model. 
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FIG. 8. Same as Fig. [T] but for Z'^r model. 
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IV. CURRENT LIMITS FROM THE TEVATRON 

The current hmits quoted from the Tevatron p, l7| and the LHC 
theoretical calculation of <j{Z') x B[Z' — )■ ^^^~)., where i = e or fi. There the decay of 
Z' includes only the SM fermions. The most updated limits from the Tevatron range from 
about 0.8 to about 1 TeV, depending on the models. Here we are going to demonstrate that 
because of the additional decay channels of the Z', including the sfermions, neutralinos, and 
charginos, the mass limits from the Tevatron are reduced by a nontrivial amount. From 
the plots shown in the last section we saw that the branching ratio into charged leptons 
decreases substantially because of the opening of the SUSY modes. Therefore, the product 
o'(Z') X B{Z' —7- £^i^) decreases. We show in Fig. [9]the limits for (a) the sequential Z' boson 
and (b) the Z^y boson for decays into SM fermions only and for decays into SM fermions and 
SUSY particles. The choice of SUSY parameters is given in Eq. (IT4|) for set (A). With these 
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FIG. 9. Exclusion limits on (a) the sequential SM Z' boson and (b) the Z'^ boson. The lines are 
the theoretical production cross section of cr(Z') x B{^Z' — t- fi^/j,^) for (i) decays into SM fermions 
only and (ii) decays into SM fermions and SUSY particles. The solid (red) line shows the most 
current upper limit from the CDF. 

choices the Z' cannot decay into squarks but can decay into slepton pairs, sneutrino pairs, 
neutralino and chargino pairs, plus the original SM fermions. Note that the Z' can also 
decay into Higgs-boson pairs such as Ah^, AH^, and H^H^ . The limit for the sequential 
Z' boson shifts from 1012 GeV down to 989 GeV (a difference of 23 GeV). The other Z' 
models behave similarly. We show the Z' mass limits for various models when the SUSY 
decay modes are open in Fig. [TOl We summarize the changes in limits in Table IIIII for the 
parameter sets (A) and (B). Compared with the case of SM fermions only, the reduction 
in the Z' mass limits ranges between 11 (5) and 32 (18) GeV for various models with the 
SUSY parameters of set (A) (set (B)). 

TABLE III. Table showing the down shift of the limits on the Z' boson mass due to the inclusion 
of the SUSY decay modes of Z' . The choices of SUSY parameters are given in set (A) and in set 
(B). 
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Z' exclusion limits for the Eg and B-L models (4.6 fb ) 
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FIG. 10. Exclusion limits on Z' bosons of various models described in Sec. II when the SUSY 
decay modes are open. The SUSY parameters chosen are listed in Eq. ()14p . The lines are the 
theoretical production cross section of a{Z') x B{Z' — > /i+/i~). The solid (red) line shows the most 
current upper limit from the CDF. 

V. PROJECTED SENSITIVITIES AT THE LHC 



Here we investigate the sensitivities at the LHC with ^/s = 7, 10, and 14 TeV with 
projected luminosities of 1, 10, and 100 fb^^, respectively, using the conventional search 
channel, pp — ?■ Z' — )■ i~^i~. Since the Z' boson is a narrow resonance state, we adopt the 
following procedures to obtain the 5cr discovery potential of each Z' model. 

1. Take a Z' boson of mass, say, 800 GeV. We pick a conservative bin-size resolution of 
order 100 GeV in the Mu distribution around TeV. Therefore, we are looking at the 
window (800 - 50, 800 + 50) GeV. 

2. We calculate the expected Drell-Yan background in this mass window with the center- 
of-mass energy and the corresponding luminosity, say A^ events. 

3. If A^ > 10 we can use Gaussian distribution, and the standard deviation is given by 
A/iV. The 5a signal corresponds to 5 x y/N events. If A^ < 10 then we should use 
Poisson statistics to determine the 5a signal events. 
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4. Check the cr{Z') x B{Z' — )> £~^i~) if it is larger or smaller than the 5a events. If 
cr(Z') X B{Z' — )■ i~^£~) is larger, it means that the discovery sensitivity can go further; 
otherwise, the discovery sensitivity cannot go further and the current Mz' is the 5a 
discovery limit, and the procedures stop here. 

5. Repeat with an increment of, say 10 GeV, for Mz' and go back to step (1). 

The 5a discovery sensitivity for LHC-7, LHC-10, and LHC-14 with luminosities 1, 10, 
and 100 fb~^, respectively, are shown in Table [IV] for the decay into SM particles only, 
and for decay into SM and SUSY particles with SUSY parameter sets (A) and (B). The 
difference in sensitivity reach between set (A) and (B) is very small here, because the mass 
of the Z' boson near the sensitivity reach is so large that the masses of the sleptons and 
gauginos involved are relatively negligible. The mass reach for the Eq Z' bosons is about 
1.3-1.5 TeV at the LHC-7 with 1 ih-\ and could be up to 2.5-2.6 TeV at the LHC-10 with 
10 fb'^, and 4.2 — 4.3 TeV at the LHC-14 with 100 fb~^ when supersymmetry is included 
in the decay of the Z'. On the other hand, the mass reach for Zgj^ and Z'q_]^ is about 
one to a few hundred GeVs better than the Z' in Eq models. We have checked that the 
number of signal events at the mass-reach limit, say for the Z'j^ model, is 5 at the LHC-7 
against a Drell-Yan background of 0.09 events, and is 4 at the LHC-10 and LHC-14 against a 
Drell-Yan background of 0.04 and 0.03 events, respectively. The probability for background 
fluctuation is < 10^^, such that it is a 5a discovery. Note that we only count one channel of 
the charged lepton (electron). If both electron and muon are counted, there will be 8 — 10 
charged-lepton pair events for discovery, which should be clean enough against negligible 
Drell-Yan background. 

VI. SUPERSYMMETRIC DECAY MODES 

Identification of supersymmetric decay modes of the Z' has it own interests, namely, to 
understand the role of Z' in the SUSY breaking. Moreover, if the Z' boson decays frequently 
into SUSY particles, we can make use of the SUSY channels to probe for the Z'. So far, in 
the models that we illustrate the branching ratio into charged leptons is not negligible, such 
that the best discovery mode is still the charged-lepton mode, which cleanly shows the peak 
in the invariant-mass distribution. Nevertheless, there exist models, e.g, Refs. [51], in which 
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TABLE IV. The 5a discovery sensitivity reach in Mz' (TeV) for various Z' models at the LHC-7, 
-10, and -14 with luminosities 1, 10, and 100 fb^^, respectively. The cases with decay into SM 
particles only, and decay into SM and SUSY particles with parameter sets (A) and (B) are shown. 





Z'n 


^x 


^; 


^; 


^scc 


Z'l 


^SM 


Z'b-l 










LHC-7 at 1 fb-^ 








SM 


1.50 


1.53 


1.46 


1.35 


1.53 


1.52 


1.83 


1.74 


SUSY [Set (A)] 


1.36 


1.49 


1.32 


1.31 


1.48 


1.46 


1.68 


1.60 


SUSY [Set (B)] 


1.37 


1.49 


1.33 


1.31 


1.48 


1.47 


1.69 


1.61 










LHC-10 at 10 fb~i 








SM 


2.65 


2.69 


2.61 


2.53 


2.67 


2.65 


3.16 


3.03 


SUSY [Set (A)] 


2.50 


2.62 


2.45 


2.45 


2.59 


2.56 


3.03 


2.82 


SUSY [Set (B)] 


2.51 


2.63 


2.46 


2.46 


2.59 


2.56 


3.04 


2.83 










LHC-14 at 100 fb-^ 








SM 


4.47 


4.51 


4.44 


4.34 


4.44 


4.41 


5.21 


5.04 


SUSY [Set (A)] 


4.25 


4.41 


4.21 


4.23 


4.33 


4.27 


5.03 


4.89 


SUSY [Set (B)] 


4.26 


4.41 


4.22 


4.23 


4.33 


4.27 


5.03 


4.89 



the charged-lepton decay mode is highly suppressed. One could also imagine that a Z' does 
not couple to fermions or sfermions but only to the Higgs sector, such that it couples solely 
to Higgs bosons and Higgsinos. In such an extreme the Z' would substantially decay into 
Higgsinos (or the physical neutralinos and charginos after mixings). In other words, the 
supersymmetric decay modes of the Z' boson could be sizable and useful for understanding 
the SUSY breaking. 

Typically, the SUSY decay modes include (i) Z' -^ it -^ £+£-x?X?, (ii) Z' -^ x?X2 ^ 
£+£~XiXi) (iii) Z' — 7- Xi^xr ~^ ^^^^^i^XiXi^ 6tc. Such leptonic modes give rise to a signa- 
ture consisting of a charged-lepton pair and large missing energies. It is clean and we can 
construct the cluster transverse mass Mt of the lepton pair and the missing energy. The 
transverse mass would indicate a broad peak structure, which is sensitive to the intermediate 
Z' boson mass. Here we only give a taste of what one can do to see the presence of the Z' 
via the supersymmetric decays. Detailed studies including various SUSY spectra and decay 
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modes, and branching ratios will be given in a future publication. 



A. Slepton-pair production 



Let us first investigate slepton-pair production in MSSM and in MSSM plus a Z' . Elec 



troweak production of e^e^ or eRe\ goes through the 7, Z, Z' exchanges 

cross section for the subprocess qq — )■ eQ,e* [a = L, R) is given by 

da 13 



B 



The differential 



dcos9 967rs 
where 



\MLa\^ + \MR^\^) {ui-mlj, (16) 



MLa - -. + T 2 , • -p 9l9 + T 2 , • -p <^<7LVe, ' UO 

s s — Tfir, -\- irrizi- z s — Tfiy, -\- luiz'i- z' 



2 I ■ -n 3r9 + - 2 I • -p 

s — m^z + ^"^zf z s — ^^z' + i^fnz'i- z' 



M«„ = ^l^^f^ + , ^, '] _ ^ gW" + l-—T^T--^Q'.nQl ' (18) 



and /3 = a/1 — 4m|^/s. Here the electric charge Qg^ = Qe, the Z charge g^°' = g^, and the 
Z' charge Q'g^ = Q'e^- The subprocess cross section is then folded with parton distribution 
functions to obtain the total cross section. The so-produced e^e^ will decay into the electron 
and positron and the lightest neutralinos under the normal hierarchy of SUSY masses, l 
Thus, the final state consists of a charged-lepton pair and a missing energy. We can construct 
the cluster transverse mass given by 

2 ^/^ 

Mt = { Xlp'r.+e- + K+.-+ l^T ] - [pTe+e- + ^T 



[yPle+e 


_+M2 


— ^ 

fe-+ ^T 
J 


2 


^Te+e- + V 


/pTe+e- 


+ K..- 


; 



(19) 

where the second equality is because pre+e- = —"^t- We show the distribution for this cluster 
transverse mass in Fig. [TTJ We have imposed a set of leptonic cuts before we construct the 
cluster transverse mass: 

PTe > 30 GeV, |r/e| < 3, ^r > 50 GeV . (20) 

The Z' models shown in Fig. [TT]are Z'j^, Z' Z' and Z'^_^. The other Z' models show similar 

features. The underneath curve is the MSSM contribution only with 7 and Z exchanges 

^ In this work, we will ignore all the Higgs exchange diagrams since they are suppressed by light quark 

masses. 
'^ Here we assume the lightest neutralino is the lightest supersymmetric particle (LSP) and the usual order 

of SUSY masses: TTig^ ~ rngj^ > m^o . 
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while the upper curve includes also the contribution from Z' . The Z' peak becomes broad 
because of the missing energies from the two neutralinos involved. Nevertheless, the sharp 
edge of the peak is sensitive to the mass difference between the Z' and the slepton masses. 
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FIG. 11. The differential cross section versus the cluster transverse mass M^ defined in Eq. (J19p for 
VP ~^ 7) Z^ Z' — )• clc^ + e_Re|j — )• e^~e~x5xi- We have applied the leptonic cuts given in Eq. (pOj) . 
Both e^^e^ and e^ejj are included. The underneath curve is without the Z' boson while the upper 
curve includes the Z' boson. The Z' models shown here are Z'^,Z'Z' and Z'q_j^. The masses 
are M^/ = 1.5 TeV, mg^^^ = 200 GeV, and m^^o = 80 GeV. 



B. Neutralino-pair production 

Next, we study the production of a neutralino pair x^^H- Assuming x\ is the LSP, then 
X2 can decay into x?^"*"^" via a virtual Z boson or a virtual slepton. The final state consists 
of a charged lepton pair and a missing energy. We can again construct the cluster transverse 
mass as in Eq. (IT^ . 

Electroweak production of X1X2 go^s through the Z and Z' exchanges. We assume that 
the squarks are much heavier such that the t-channel squark exchanges are suppressed. The 
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differential cross section for the subprocess q{pi)q{p2) — > Xi(^i)X2(^2) is given by 
da P' 



dcosO 96vrnKs-^)K'-^) (|M..(S)P + |M,,(5)P) 
+(m|o - i) (m|o - i) {\Mlr{s)\^ + \MnLm') 
+2 m^o m.^0 s (Re [Mll{s) M^^(s)] 



+Re[MnRis)MUs)])Y (21) 

where 

Ma,is) = -. /] . ^ Ol^gj + -^ ^. —C';-Q'^, , (22) 

a, (3 = L,R and 

/3'=|[l-(m|o + m|o)/s] -(2m^om^o/s)2| . (23) 



Here the chiral couphngs 0-^2 of the Z boson and the chiral couphngs C-^^2 of ^^e Z' boson 
to the neutrahnos are, respectively, given by 

o;? = -o:t , (24) 

and 

c;;^ = -c';^'* , (25) 

where N is the mixing matrix of the neutrahnos defined in the appendix. Numerically, 
with the choice of SUSY parameters of set (A), we obtain the masses m^o = 74.8 GeV and 
m^o = 130.6 GeV, and the mixing parameters Ni^ = 0.52, A'^23 = —0.41, A'^14 = —0.33, and 
N24 = 0.38. 

We apply the same set of leptonic cuts as in Eq. (!20l) and construct the cluster transverse 
mass. We show the cluster transverse-mass spectrum in Fig. [T2j The Z' models shown in 
Fig. [12] are Z^, Z'^, Z'^, and Z^^^. The underneath curve is the MSSM contribution only 
with the Z exchange while the upper curve includes also the contribution from Z'. The 
Z' peak becomes broad because of the missing energies from the two neutrahnos involved. 
Nevertheless, the edge of the peak is sensitive to the mass difference between the Z' and the 
neutrahnos. 
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FIG. 12. The differential cross section versus tlie cluster transverse mass Mt defined in Eq. (J19p for 
pp -^ Z, Z' ^ X1X2 followed by the leptonic decay of X2 ~^ Xie~'"e~. We have applied the leptonic 
cuts given in Eq. (j20p and assumed the branching ratio i?(x2 ~^ XiC''"e~) = 0.1. The underneath 
curve is without the Z' boson while the upper curve includes the Z' boson. The Z' models shown 
here are Z^, Z!^, Zj, and Z!,^^. The masses are Mz' = 1.5 TeV, m^o = 74.8 GeV, and m^o = 130.6 
GeV. 

C. Chargino-pair production 

Lastly, we have electroweak production of xtXi that goes through the 7, Z, and Z' 
exchanges. The differential cross section for the subprocess q{pi)q{p2) — > Xi^(^i)xr(^2) is 
given by 



da 



/3 



dcosO 



+Re[MnRirs)Ml^{s)])Y 



(26) 
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where 



M„,(,) .f!«^ . -4f^ . . f'f'^ , (27) 

s s — m^ + imzi- z s — m^, + zmz'i z' 



a, 13 = L, R, and (3 = . h- — Arni^/s. Here the chiral couphngs 0{^ of the Z boson and the 
couphngs C^ of the Z' boson to charginos are, respectively, 

On = |V^iir + ^1^121' -a:w, 

On = |f/ii|' + ^|f/i2r-a;„, (28) 

and 

On = Q//„l^i2| , 

C[^, = -Q'hPi2\\ (29) 

where U and V are the mixing matrices of the charginos. Numerically, with the choice in 
set (A) for SUSY parameters we have ^12 = 0.75, Vu = -0.67, Uu = 0.89, Uu = -0.46 
and m-+ = 109.5 GeV. 

We only calculate the production of xtXi ^ because the second chargino is about twice as 
heavy as the first one. Each of the charginos decays via a virtual W, z/, or i into a charged 
lepton, a neutrino, and the lightest neutralino (if going through the virtual W, light quarks 
are also possible). Therefore, there will two charged leptons plus missing energies in the 
final state. Just as the same as the case of slepton-pair or neutralino-pair production, we 
reconstruct the cluster transverse mass as in Eq. (fT9|) . We show the distribution of cluster 
transverse mass in Fig. [13] for Z'j^, Z'^, and Z'^^^. It is easy to see the bump due to the 
presence of the Z' boson, though the bump is not as discernible as the previous two cases 
of slepton-pair and neutralino-pair production. 

VII. CONCLUSIONS 

In this work, we have studied the possible supersymmetric decay modes of an additional 
neutral gauge boson Z', which is currently limited to be at least 1 TeV. Grand unified 
theories have predicted one or more such Z' bosons along the path through which the 
GUT symmetry is broken down to the electroweak symmetry. When supersymmetry is 
included in the theory, such a Z' boson can decay not only into the SM particles but 
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FIG. 13. The differential cross section versus the cluster transverse mass Mt defined in Eq. (J19p 
for pp — ;• xtxi followed by the leptonic decay of Xi -^ Xi^^^e- We have applied the leptonic cuts 
given in Eq. (j20p and assumed the branching ratio B[xt ~^ Xi^^^e) x B{xi -^ x!i^~^e) = 0.1. 
The underneath curve is without the Z' boson while the upper curve includes the Z' boson. The 



Z' models shown here are Z'^, Z' and ^scc- '^^^ masses are Mz' = 1.5 TeV and m- 
GeV. 
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also the supersymmetric partners. We have used E^, U{1)b-l, and the sequential models 
to illustrate how the decays of the Z' are affected. In particular, the golden search mode- 
charged leptons- for the Z' will have a smaller branching ratio as the supersymmetric modes 
open. We have shown that the current limits obtained at the Tevatron and the LHC will be 
reduced by a noticeable amount, of order 20 GeV. We have also estimated the 5a discovery 
sensitivities of the Z' at the LHC, including the effect of supersymmetric decay modes. 
Finally, we demonstrated that even though the Z' decays into supersymmetric particles, 
giving rise to missing energies, one can still reconstruct the cluster transverse mass (using 
the observable charged leptons) to identify the existence of the Z'. We believe further studies 
along this direction is worthwhile since it can help us to fully understand the role of the Z' 
in the supersymmetry-breaking and the symmetry breaking pattern. 
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Appendix A: Feynman Rules 

As explained in Sec. IIB, the extra term hgSHuHd allowed by the U{iy symmetry will 
generate the effective yU. term when the singlet scalar field S develops a VEV. With the 
assumption that the singlet scalar and singlino fields, Z'-ino, and other exotic fermions are 
very heavy, the low-energy particle content includes the Z' boson and those of MSSM. The 
effective superpotential Wes involving the matter and Higgs superfields is the same as the 
MSSM's, and given by 

W^eff = e,t [y-jQ^jHtU^ - yiQ^jH'.DI - yl^L]H^E^ + ^.H^H^] , (Al) 

where e^ = —621 = ^^hj are family indices, and y" and y'^ represent the Yukawa matrices for 
the up-type and down-type quarks respectively. Here Q, L, f/^, D'^, E^, Hu, and Ha denote the 
superfields for the quark doublet, lepton doublet, up-type quark singlet, down- type quark 
singlet, lepton singlet, up- type Higgs doublet, and down-type Higgs doublet respectively. 
The scalar interactions are obtained by calculating the F- and D-terms of the superpotential, 
and by including the following soft-SUSY-breaking terms 

£soft = -MiBB - M2WW - M^gg 

- M^Q^Q - Ml^U'^^W - Ml^D^^D"" - MIDl - M^E'^^E'' (A2) 

Am = ea6 [ylAlQ^^H'Jj^ - yf^A^^Q]H',D1 - y\^Af^^H',E<l + ^^BH:H'^ , (A3) 

where the £soft represents the soft mass terms for the gauginos l| and sfermions, and £trii 
represents the trilinear A terms. 

The gauge interactions for the fermionic and scalar components, denoted generically by 
ip and respectively, of each superfield mentioned above are given by 

C = i, i^D^ ^ + {D^y {D^<f)) , (A4) 

^ There should also be a soft gaugino mass term for the Z'-ino. However, we will assume it is heavy and 
let it decouple from the low energy spectrum in the present work. 
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where the covariant derivative is defined as usual 



D^ = d^ + leQA^ + '^-4(t+W^+ + T-W-) + ig^{T^L - x^Q)Z^ + ig^Z'^Q . (A5) 

Here g\ = g/ cos9„ with g the SU{2)l gauge couphng and 6^ the weak mixing angle, g2 is 
the gauge coupling for the extra f/(l), r^ and T^^ are the ladder operators and the third 
component of the SU{2)l generators, Q and Q' are the charges of the two f/(l)s in unit 
of the electromagnetic charge e and g2 respectively, and finally, x„ = sin^6'w. Since we 
only consider supersymmetric f/(l) symmetry for the Z' boson and due to the Majorana 
nature of the Z'-ino, there is no coupling between the Z'-ino and the Z' boson. Simply from 
Eqs. (lA4p and (lASp we obtain the interactions of Z' with fermions, sfermions, neutral and 
charged Higgsinos (which become the physical neutralinos and charginos after mixing effects 
are taken into account), and the Higgs bosons. 

The rotation of neutral bino, wino, and Higgsinos into the physical neutralinos is given 
by 



\xlj 



N 






(A6) 



where A^ is an orthogonal matrix. The rotation of the charged wino and Higgsino into the 
physical charginos is via bi-unitary transformation 



U 



(A7) 



xl \ ^, I W^ \ I xi 

X2 I \ K I \X2 / \ -d 

where U and V are unitary matrices. The mixing angles involved in physical Higgs bosons 
{h^, H^, H^, A^) can be read off from the following decompositions of the Higgs boson fields 
Hu and Ha. 



H,, 



H, 




H+cos(3 + G+sin(3 
^ + ^(/i° cos a + i/o sin a) + z(AO cos /3 + G° sin /3) 

^^ + 75 (-^° sin a + H^ cos a) + t{A^ sin P - G^ cos /3) 



(A8) 



(A9) 



H sin/3 — G cos/3 

where the angle a is the mixing of the neutral CP-even Higgs bosons /i" and iiT", G°'^ are 
the Goldstone bosons, and v = 246 GeV. 
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Lastly, there are Yukawa-type interactions between the gauginos and the scalar and 
fermionic ip components of the matter superfield of the following form 

C = -y/2gy*T''X''iP + h.c. (AlO) 

where a is the group index of the f/(l)y, SU{2)l, SU{3)c, or the extra U{l)z'- The inter- 
actions involving the Z'-ino are 

C^, = -v^(^2 Q^Z'Q'qQ + U'^Z'Q'uM" + D'^^Z'Q'^.D'' + h.c. 
+L^Z'Q'lL + E'^^Z'Q'e.E^ + h.c. 
+HlZ'Q'H^H^ + H\Z'Q'Ha + \i.c. . (All) 



Note that the Z'-ino will mix with the H^ and Hd Higgsinos when the /7° and H^ take on 
vacuum expectation values. Thus, we will have a 5 x 5 neutralino mass matrix. However, 
we decouple the Z'-ino in this work by setting the Z'-ino mass heavy. We will come back to 
this in later work. 

In the following we list the Feynman diagrams and the corresponding Feynman rules 
involving the Z' boson and the MSSM particles. For each model, one should use the corre- 
sponding coupling strength and chiral charges. The coupling strength g2 and the Z' charges 
are given in Tables HI and HTl for the Eq and U{1)b-l models respectively. For the sequential 
Z' model, replace the coupling strength g2 by gi while the corresponding chiral charges are 
given by the SM values: Q'^^Z'^^j) = T^i — x^Qi. The chiral couplings of the charginos and 
neutralinos with the Z' boson are, respectively, given by 

cl' = -q'h,n,sn;, - q'^^n,,n;, , C:^ = -cll^* . (A12) 

These coupling coefficients are the same for the three Z' models that we have studied in this 
work, as long as we use the corresponding Q' charges for the two Higgs doublet fields. 
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